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The most general method for encoding quantum information is not to encode the information 
into a subspace of a Hilbert space, but to encode information into a subsystem of a Hilbert space. 
Recently this notion has led to a more general notion of quantum error correction known as operator 
quantum error correction. In standard quantum error correcting codes, one requires the ability to 
apply a procedure which exactly reverses on the error correcting subspace any correctable error. In 
contrast, for operator error correcting subsystems, the correction procedure need not undo the error 
which has occurred, but instead one must perform correction only modulo the subsystem structure. 
This does not lead to codes which differ from subspace codes, but does lead to recovery routines 
which explicitly make use of the subsystem structure. Here we present two examples of such operator 
error correcting subsystems. These examples are motivated by simple spatially local Hamiltonians 
on square and cubic lattices. In three dimensions we provide evidence, in the form a simple mean 
field theory, that our Hamiltonian gives rise to a system which is self-correcting. Such a system 
will be a natural high-temperature quantum memory, robust to noise without external intervening 
quantum error correction procedures. 

PACS numbers: 03.67.Lx,03.67.a,75.10.Pq 



I. INTRODUCTION 



In the early days of quantum computation, the analog nature of quantum information and quantum transforms, 
as well as the effect of noise processes on quantum systems, were thought to pose severe obstacles^, Q towards the 
experimental realization of the exponential speedups promised by quantum computers over classical computers H, EL 
Soon, however, a remarkable theory of fault-tolerant quantum computation]?!. Isl lij. ITol ITU Il2l ll3L FLU ITU Il6l Il7| 
emerged which dealt with these problems and showed that quantum computers are indeed more similar to probabilistic 
classical computers than to analog devices. Analog computers have a computational power which is dependent on a 
lack of noise and on exponential precision, whereas probabilistic classical computers can be error corrected and made 
effectively digital even in the presence of noise and non-exponential precision. The theory of fault-tolerant quantum 
computation establishes that quantum computers are truly digital devices deserving of the moniker computer. An 
essential idea in the development of the theory of fault-tolerant quantum computation was the notion that quantum 
information could be encoded into subspaces0, 01 |jj 0] (quantum error correcting codes) and thereafter protected 
from degradation via active procedures of detection and correction of errors. Encoding quantum information into 
subspaces, however, is not the most general method of encoding quantum information into a quantum system. The 
most general notion for encoding quantum information is to encode the information into a subsystem of t he q uantum 
system^, 0, This has been perha ps b est exploited in the theory of noiseless subsystems pH I20L Ell I22L l23t |24 | 
and dynamic recoupling schemes |22ll25ll26i |. Recently a very general notion of quantum error correction has appeared 
under the moniker of "operator quantum error correction." [27l |2S| In this work the possibility of encoding into 
subsystems for active error correction is explicitly examined. While it was found [27ll28l| that the notion of a encoding 
into a subsystem does not lead to new codes (all subsystem codes could be thought of as arising from subspace codes), 
encoding into a subsystem does lead to different recovery procedures for quantum information which has been encoded 
into a subsystem. Hence operator quantum error correcting codes, while not offering the hope of more general codes, 
do offer the possibility of new quantum error correcting routines, and in particular to the possibility of codes which 
might help improve the threshold for fault-tolerant quantum computation due to the lessened complexity of the error 
correcting routine. 

In this paper we present two examples of operator quantum error correcting codes which use subsystem encodings. 
The codes we present have the interesting property that the recovery routine does not restore information encoded 
into a subspace, but recovers the information encoded into a subsystem. Using the [n, d, k] labelling a quantum 
error correcting codes, where n is the number of qubits used in the code, d is the distance of the code, and k is the 
number of encoded qubits for the code, our codes are [n 2 ,n, 1] and [n 3 ,n, 1] quantum error correcting codes. The 
subsystem structure of our codes is explicitly exploited in the recovery routine for the code, and because of this they 
are substantially simpler than any subspace code derived from these codes. 

While the two codes we present are interesting in there own right, there is a further motivation for these codes above 
and beyond their exploitation of the subsystem structure in the recovery routine. The two operator quantum error 
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correcting subsystems we present are motivated by two interesting Hamiltonians defined on two and three-dimensional 
square and cubic lattices of qubits with certain anisotropic spin-spin interactions |2 1| . The three-dimensional version 
of this system is particularly intriguing since it offers the possibility of being a self-correcting quantum memory. In 
a self-correcting quantum memory, quantum error correction is enacted not by the external control of a complicated 
quantum error correction scheme, but instead by the natural physics of the device. Such a quantum memory offers 
the possibility for removing the need for a quantum microarchitecture to perform quantum error correction and could 
therefore profoundly speed up the process of building a quantum computer. In this paper we present evidence, in the 
form of a simple mean field argument, that the three-dimensional system we consider is a self-correcting quantum 
memory. We also show that the operator error correcting subsystem structure of this code is an important component 
to not only the self-correcting properties of this system, but also to encoding and decoding information in this system. 

The organization of the paper is as follows. In Section [H] we review the notion of encoding information into a 
subsystem and discuss the various ways in which this has been applied to noiseless subsystems and dynamic recoupling 
methods for protecting quantum information. Next, in Scction lTTTl we discuss how operator error correcting subsystems 
work and how they differ from standard quantum error correcting codes. Our first example of a operator quantum 
error correcting subsystem is presented in Section Hvl where we introduce an example on a square lattice. The second, 
and more interesting, example of a operator quantum error correcting subsystem is given in Section where we 
discuss an example on a cubic lattice. In Section IVT1 we introduction the notion of a self-correcting quantum memory 
and present arguments that a particular Hamiltonian on a cubic lattice related to our cubic lattice subsystem is self- 
correcting. We conclude in Section IVIII with a discussion of open problems and the prospects for operator quantum 
error correcting subsystems and self-correcting quantum memories. 



II. SUBSYSTEM ENCODING 



Consider two qubits. The Hilbert space of these qubits is given by C 2 ® C 2 . Pick some fiducial basis for each 
qubit labelled by |0) and One way to encode a single qubit of information into these two qubits is to encode the 
information into a subspace of the joint system. For example, we can define the logical basis states |0l) = "^(|01) — 

10}) and = 1 11) such that a single qubit of information can be encoded asa|0i)-|-/3|lL) = ^=(|01) — 110)) +/3|11). 
This is an example of idea of encoding quantum information into a subspace, in this case the subspace spanned by 
|0l) and |1l). But another way to encode a single qubit of information is to encode this information into one of the 
two qubits. In particular if we prepare the state \ip) (g) (a|0) +/3|1)) for an arbitrary single qubit state then we have 
also encoded a single qubit of information in our system. This time, however, we have encoded in the information 
into a subsystem of the system. It is important to note that the subsystem encoding works for an arbitrary state 
\ip). If we fix \tp) to some known state, then we are again encoding into a subspace. We reserve the nomenclature of 
"encoding into a subsystem" to times in which is arbitrary. 
More generally, if we have some Hilbert space Ti., then a subsystem C is a Hilbert space arising from Ji as[9l Il8| 

H = {C®V)®£. (1) 

Here we have taken our Hilbert space and partitioned it into two subspaces, £ and a subspace perpendicular to £. On 
this perpendicular subspaces, we have introduced a tensor product structure, C®T>. We can then encode information 
into the first subsystem C. This can be achieved by preparing the quantum information we wish to encode pc into 
the first subsystem, C along with any arbitrary state po into the second subsystem T>: 

P=(Pc®Pd)®0- (2) 

The fact that quantum information can most generally be encoded into a subsystem was an essential insight used 
in the construction of noiseless (decoherence-free) subsystems 0, 113, El HH HH 13 ■ Suppose we have a system 
with Hilbert space Tis and an environment with Hilbert space 'He- The coupling between these two systems will be 
described by an interaction Hamiltonian H int which acts on the tensor product of these two spaces Hs ®'Hb- The 
idea of a noiseless subsystem is that it is often the case that there is a symmetry of the system-environment interaction 
such that the action of the interaction Hamiltonian factors with respect to some subsystem structure on the system's 
Hilbert space, 

Him =^[{Id®D a )@E a )®B a , (3) 

a 

where 1^ is the d-dimensional identity operator acting on the subsystem code space C, D a acts on the subsystem V, 
E a acts on the orthogonal subspace £, and B a operates on the environment Hilbert space He- When our interaction 
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Hamiltonian possesses a symmetry leading to such a structure, then, if we encode quantum information into C, this 
information will not be effected by the system-environment coupling. Thus information encoded in such a subsystem 
will be protected from the effect of decoherence and hence exists in a noiseless subsystem. Noiseless subsystems were 
a generalization of decoherence- free subspaces[2|| [3(3, this latter idea occurring when the subsystem structure is not 
exploited, T> = C, and then encoding quantum information is simply enco ding quantum information into a subspace. 
Subsystems have also been used in dynamic recoupling techniques p3. l2fi| where symmetries are produced by an 
active symmetrization of the system's component of the system-environment evolution. 



III. OPERATOR QUANTUM ERROR CORRECTING SUBSYSTEMS 

Here we examine the implications of encoding information into a subsystem for quantum error correcting 
protocols[2?l |2S| . Suppose that we encode quantum information into a subsystem C of some quantum system with 
full Hilbert space H — (C ® X>) © £. Now suppose some quantum operation (corresponding to an error) occurs on 
our system. Following the standard quantum error correcting paradigm, we then apply a recovery procedure to the 
system. When D = C, i.e. when we are encoding into a quantum error correcting subspace, then the quantum error 
correcting criteria is simply that the effect of the error process followed by the recovery operation should act as identity 
on this subspace. If we encode information into a subsystem, however, this criteria is changed to only requiring that 
the recovery operation should act as identity on the subsystem C. In particular we do not care if the effect of an error 
followed by our recovery procedure enacts some nontrivial procedure on the T> subsystem. In fact our error correcting 
procedure may induce some nontrivial action on the T> subsystem in the process of restoring information encoded in 
the C subsystem. 

How does the above observation modify the basic theory of quantum error correcting codes? In standard quantum 
error correction, we encode into some error correcting subspace with basis |z). The necessary and sufficient condition 
for there to be a procedure under which quantum information can be restored under a given set of errors E a is given 
by 

(i\ElE b \j) = 6 itj c a , b . (4) 

For the case of encoding into a subsystem this necessary and sufficient condition is modified as follows. Let \i) ® \k) 
denote a basis for the subspace C ® V. Then Kribs et aZ. |27ll2g| showed a necessary condition pAllCl] for the quantum 
error correcting is given by 

((t| ® (k\) ElE b (|i) ® |/)) = 5 itj m a , b , k ,l- (5) 

That this condition is also sufficient has recently also been shown by Nielsen and Poulin|3l]|. As noted in |27l |28[. a 
code constructed from the subsystem operator quantum error correcting criteria can always be used to construct a 
subspace code which satisfies the subspace criteria Eq. @. We note, however, that while this implies that the notion 
of using subsystems for quantum error correction does not lead to new quantum error correcting codes above and 
beyond subspace encodings, the codes constructed which exploit the subsystem structure have error recovery routines 
which are distinct from those which arise when encoding into a subspace. In particular, when one encodes into a 
subsystem, the recovery routine does not need to fix errors which occur on other subsystems. Below we will present 
examples of subsystem encodings in which the subsystem structure of the encoding is essential not for the existence 
of the quantum error correcting properties, but it essential for the simple recovery routine we present. 



IV. TWO-DIMENSIONAL OPERATOR QUANTUM ERROR CORRECTING SUBSYSTEM 

Here we construct an operator quantum error correcting subsystem for a code which lives on a two-dimensional 
square lattice. This code makes explicit use of the subsystem structure in its error recovery procedure. A familiarity 
with the stabilizer formalism for quantum error correcting codes is assumed (see [32113 for overviews.) 



A. Preliminary Definitions 



Consider a square lattice of size n x n with qubits located at the vertices of this lattice. Let C\j denote the operator 
O acting on the qubit located at the ith row and jth column of this lattice tensored with identity on all other qubits. 
Recall that the Pauli operators on a single qubit are 



X 



1 

1 



,iY = 



-1 

1 



and Z = 



1 

-1 



(6) 
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It is convenient to use a compact notation to denote Pauli operator on our n qubits by using two n bit strings, 

p( 0l 6) = n x id z it = n i ^ if fl « = ° and ''■.< 1 • w 

»,j=i »,i=i [ — if Oij = 1 an d hj = 1 

where a, 6 6 are n by n matrices of bits. Together with a phase, i 9 *, <p 6 Z4, a generic element of the Pauli group 
on our n 2 qubits is given by i^P(a, b). We will often refer to a Pauli operator as being made up of X and Z operators, 
noting that when both appear, the actual Pauli operator is the iY operator. 

We begin by defining three sets of operators which are essential to understanding the subsystem structure of our 
qubits. Each of these sets will be made up of Pauli operators. 

The first set of Pauli operators which will concern us, T, is made up of Pauli operators which have an even number 
of Xij operators in each column and an even number of Zij operators in each row: 

{n n j 

(-1)^(0,6)10 6 23,00^=0, and06 iJ = ol, (8) 
i=l ]=1 J 

where denotes the binary exclusive-or operation (we use it interchangeably with the direct sum operation with 
context distinguishing the two uses.) Note that these operators form a group under multiplication. This group can 
be generated by nearest neighbor operators on our cubic lattice 

T = (XijXi + ij, Zj^Zj y i + i,\/i 6 Z n _x,j 6 Z n ) . (9) 

Examples of elements of the group T are diagrammed in Fig. ^ 
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FIG. 1: Above we have represented elements of the T . Each set of operators enclosed in a rectangle represents a Pauli operator 
acting on the particular qubits tensored with the identity on all other qubits. Each of the operators enclosed in the dotted 
rectangles are elements of T . 

The second set of Pauli operators we will be interested in is a subset of T, which we denote S. S consists of Pauli 
operators which are made up of an even number of rows consisting entirely of X operators and an even number of 
columns consisting entirely of Z operators, 



{n I n \ n / n 

i=l \j=l J 0=1 \i=l 



hjj = j (10) 

where /\ is the binary and operation. S is also a group. In fact it is an Abelian subgroup of T. Further all of the 
elements of S commute not just with each other, but with all of the elements of T. It can be generated by nearest 
row and column operators, 

In n \ 

S = ( 1 { -W<-V ; ' II Z, ,Z, j 1 . : / e Z„_i ) . (11) 

\i=l i=l I 



These generators will be particularly important for us, so we will denote them by 

n n 

Sf = TT Xi t jXi+t t j, and = TT ZijZi^+i. 
j=i i=i 



(12) 



S is a stabilizer group familiar from the standard theory of quantum error correcting codes. An example of an element 
in S is given in Fig. [2] 



X X X 



■x—:^-x 



x—x—x 



FIG. 2: A nontrivial element of the group <S. This element has an even number of columns which are entirely X operators 
multiplied by an even number of rows which are entirely Y operators. Notice how the Y elements appear where both of these 
conditions are met. 

The final set of operators which we will consider, £, is similar to S except that the evenness condition becomes an 
oddness condition, 



n I n \ n / n \ 

C={(-l)t>P(a,b)\cP£Z 2 ,Q}lf\a i j} =1,0 =1 

j=l y=l / J=l Vi=l / 



(13) 



This set does not by itself form a group, but together with S it does form a group. This combined group is not 
Abelian. C has the property that all of its elements commute with those of T and S. A nontrivial element of C is 
given in Fig. [3] 



X X X - 



- X X X X X 



FIG. 3: A nontrivial element of the set £.. This element has an odd number of columns which are entirely X operators 
multiplied by an oddnumber of rows which are entirely Y operators. It represents an encoded Y operator on the encoded qubit 
as described in Sec llVfJl 
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B. Subsystem Structure 



We will now elucidate how T, S, and C are related to a subsystem structure on our n qubits. Let H = (C 2 ) n 
denote the Hilbert space of our n 2 qubits. 

We first note that since S consists of a set of mutually commuting observables, we can use these observables to 
label subspaces of TL. In particular we can label these subspaces by the 2(n — 1) ±l-valued eigenvalues of the Sf- 
and Sf operators. Let us denote these eigenvalues by sf and sf respectively and the length n — 1 string of these ±1 
eigenvalues by s x and s z . We can thus decompose TC into subspaces as 

H= H s x tS z = H s x >s z. (14) 

s x ,...,s x _ 1 ,sf ,...,s z _ 1 =±l s x ,s z 

By standard arguments in the stabilizer formalism, each of the H s x sz subspaces is of dimension d = 2" _2 (™ _:1 . J u st 
for completeness, we note that the operators Sf and Sf act under this decomposition as 

Sf = sf 7 2 „2_ 2( „_ 1) , 

s x ,s z 

sf = »fV-w ( 15 ) 

Now examine the two groups T and the group generated by elements of C and S. Both of these groups are non- 
Abelian. All of the elements of T and C commute with elements of S. Further, all of the elements of T and C 
commute with each other. This implies, via Schur's lemma [34l l35| . that C and T must be represented on TL s x s z by 
a subsystem action. In particular the full Hilbert space splits as 

n = Hj XjsZ ®Hfx tsZ , (16) 

s x ,s z 

such that operators from T £ T act on the first tensor product 

T = 7> !sZ ®/ 2 ,vreT, (17) 

s x ,s z 

and the operators from L G C act on the second tensor product 

L = / 2( „_i,2 ® L s x , s z ,Vi e L. (18) 

s x ,s z 

Here we have assigned dimensions and 2 to these tensor product spaces. To see why these dimensionalities 

arise we appeal to the stabilizer formalism. We note that modulo the stabilizer structure of S, C is a single encoded 
qubit. Similarly if one examines the following set of (n — l) 2 operators from T, 

j 

— Xi j — X^X n _x^, (19) 

fe=l 

where i S Z n _i, j 6 Z n _i, one finds that modulo the stabilizer they are equivalent to (n — l) 2 encoded Pauli operators. 

The subsystem code we now propose encodes a single qubit into the Hilbert space TCf x s z with sf=sf — +1, Vi £ 
Z n _i (choices with other ±1 choices form an equivalent code in the same way that stabilizer codes can be chosen 
for different stabilizer generator eigenvalues.) The code we propose thus encodes one qubit of quantum information 
into a subsystem of the n 2 "bare" qubits. We stress that the encoding we perform is truly a subsystem encoding: 
we do not care what the state of the Ti.Jx s z subsystem is. For simplicity it may be possible to begin by encoding 
into a subspace which includes our particular subsystem (i.e. by fixing the state on T~tJ x s z) but this encoding is not 
necessary and indeed, after our recover routine for the information encoded into T~ifx s z , we will not know the state of 
the Tijx s z subsystem. We will denote the Hilbert space 7if x s z with sf = sf = +l,Vi E I* n -i by Wf x=s z = { +1 } n -i ■ 
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C. Subsystem Error Correcting Procedure 

If we encode quantum information into the subsystem 7i^x— a z—r + iy n -i , then what sort of error correcting prop- 
erties does this encoding result in? We will see that the S operators can be used to perform an error correcting 
procedure which restores the information on Hfx =a z = i+i\n-i> but which often acts nontrivially on the subsystem 
"HgX—gZ— i + iyn-i ■ This exploitation of the subsystem structure in the correction procedure is what distinguishes our 
subsystem operator quantum error correcting code from standard subspace quantum error correcting codes. 

Suppose that a Pauli error P(a,b) occurs on our system. For a Pauli operator P(a, 6), define the following error 
strings: 

n n 

ej(«)=©«U. fi(b)=^h,j. (20) 

i=l j=l 

Notice that if ej = fi = Q,Vi,j, then this implies that P{a,b) is in the set T. Further note that in this case, the effect 
of P(a, b) is to only act on the T~iJ x s z subsystems, i.e. P(a, b) is block diagonal under our subsystem decomposition, 
Eq. JTBJl, acting as 

P(a,b)= E s x >s z{a,b)® I 2 . (21) 

s x ,s z 

where E s x s z (a, b) is a nontrivial operator depending on the subspace labels s x , s z and the type of Pauli error (a, b). 
Therefore errors of this form (ej = fi = 0) do not cause errors on our information encoded in 'Hf x=s z = ^ +1 j n -i ■ With 
respect to the errors of this form, the information is encoded into a noiseless subsystem[18|. 

Returning now to the case of a general P(a, 6), from the above argument we see that if we can apply a Pauli operator 
Q(c, d) such that Q(c, d)P(a, b) is a new error, call it R(a', b'), which has error strings e'j (a!) — fi(b') — 0, Vi, j, then we 
will have a procedure for fixing the error P(a, b), modulo the subsystem structure of our encoded quantum information. 
In other words, our error correcting procedure need not result in producing the identity action on the subspace labelled 
by sf = sf = +l,Vi G Z n _i, but need only produce the identity action on the subsystem 'Hfx =s z = ^ +1 y n -i ■ We can 
perform just such a procedure by using the elements of S as a syndrome for which errors of small enough size can be 
corrected. 

To see how this works, suppose P(a,b) occurs on our system. Then note that measuring Sf is equivalent to 
determining 

n 

0(6ij © bi+u) = Mb) © f l+1 (b), (22) 

3=1 

and similarly measuring Sj is equivalent to determining 

n 

0(ai,j © <H,j+i) = ej(o) © e i+ i(a). (23) 

i=l 

Note that all 2 [n — 1) of these measurements can be performed simultaneously since the elements of S all commute 
with each other. We wish to use these measurement outcomes to restore the system to e 3 (a) = fi(b) — (if possible.) 

To see how to do this, treat the fi(b) as a n bit codeword for a simple redundancy code (i.e. the two codewords 
are fi(b) — 0,Vi and /,(&) = l,Vi). A similar procedure will hold for the ej(6). Measuring the n — 1 operators Sf is 
equivalent to measuring the syndrome of our redundancy code. In particular we can use this syndrome to apply an 
error correcting procedure for the fi(b) bit strings. The result of this correction procedure is to restore the system 
to either the codeword fi(b) = 0,Vi or the codeword fi(b) = l,Vi. The former corresponds to an error correction 
procedure which can succeed (given that an equivalent procedure for the ej(a) bit strings also succeeds), whereas 
the latter procedure is one where the error correction procedure will fail. Notice that our error correcting procedure, 
when it succeeds, is only guaranteed to restore the system to fi(b) — and ej(b) = 0, and thus the full effect of the 
procedure may be to apply some nontrivial operator to the Ti.J x s z subsystems. 

Let us be more detailed in describing the error correcting procedure for the fi(b) code words. Let sf be the result of 
our measurements of the Sf operators. Give the sf we can construct two possible bit strings f- and -if I (-> denotes 
the negation operation) consistent with these measurements. Let H(f') and H(—if) denote the Hamming weight of 
these bit strings (i.e. H(f') is the number of Is in the n bits f[) and define /" to be the bit string /' or -if with the 



8 



smallest Hamming weight. We now apply an operation consisting only of Zij operators. In particular we apply the 
operator 

n s 

Qi{f') = Y[Z^\ (24) 

i=l 

for any fixed column index jo. The operator Q±(f")P(a,b) is then seen to be of one of two forms: either this new 
operator has the Z error string equal to all zeros or all ones. In the first case we have successfully restored the system 
to the all fiib) = codeword, whereas for the second case, we have failed. How many Z errors can be corrected in 
this fashion? If P(a, b) consisted of Z errors b with an error string fi(b) with a Hamming weight of this string H(f) 
which is less than or equal to L'^y^Jj then the correction procedure will succeed. Thus the code we have constructed 
is a [n 2 ,n, 1] code: it encodes a single qubit into n 2 qubits and has a distance n. 

Above we have focused on the case of Pauli Z errors. Clearly an analogous argument holds for Pauli X errors 
(with the role of the rows and columns reversed.) Further, Pauli Y errors are taken care of by the combined action 
of these two procedures. By the standard arguments of digitizing errors in quantum error correcting codes, we have 
thus shown how our operator quantum error correcting subsystem code can correct up to L'^Y^J arbitrary single qubit 
errors. 



D. Logical Operators 

We comment here on the logical operators (operators which act on the encoded subsystem) for this code. From our 
analysis of the subsystem structure, it is clear that elements of £ act on the subsystem. Thus, for instance, the effect 
of a row of Pauli X operators is to enact an encoded Pauli X operation on the coded subsystem. We can choose a 
labelling of the subsystem such that 

n 

X = I[Xi,j= (25) 

j = l s x ,s z 

while in the same basis the effect of a column of Pauli Z operators is to enact and an encoded Pauli Z on the coded 
subsystem, 

n 

2 = n z m = v-i> 2 ® z - ( 2e ) 

i=l s x ,s z 

These two operators can then be used to enact any Pauli operator on the encoded quantum information. Notice 
that other elements of C also act as encoded Pauli operators on W^x =s z = / + un-i (but act with differing signs on the 

other s x ,s z labelled subspaces.) An important property of these logical operators is that they can be enacted by 
performing single qubit operators and no coupling between the different qubits is needed. This is important because 
it will allow us to assume an independent error model for error which occur when we imprecisely implement these 
gates on our encoded quantum information. Not only can the above construction be used to implement the Pauli 
operators on our subsystem code, it can also be used to measure the Pauli operators on our subsystem code. 

Another easily implementable operation on our code is a logical controlled-not. Suppose we take two identically 
sized two-dimensional codes and stack them on top of each other. Then the application of a transverse controlled-not 
operator between all n 2 of these two systems will enact a logical controlled-not between the two encoded qubits. To 
see this note that if we treat the elements of the set 5 as a stabilizer code, then these transverse operators preserve 
the combined stabilizer S x S and that the action of the n 2 controlled-not gates do not mix the C x C and T x T 
operators. 

Gottesman0> Ei3 nas shown that given the ability to measure and apply the encoded Pauli operators along with 
the ability to perform a controlled NOT on a stabilizer code, one can perform any encoded operation which is in the 
normalizer of the Pauli group (i.e. the gate set relevant to the Gottesman-Knill theoremjs^). 

We have seen how to implement encoded Pauli operators and the controlled-NOT on the information encoded into 
our subsystem. These operations do not allow for universal quantum computation, so an important open question for 
our subsystem code is to find an easily implementable method for completing this gateset to a universal set of gates. 
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E. Hamiltonian Model of the two-dimensional Subsystem Code 

An interesting offshoot the above two-dimensional operator quantum error correcting subsystem is the analysis of 
a particularly simple Hamiltonian whose ground state has a degeneracy which corresponds to the subsystem code. 
We introduce this Hamiltonian here in order to make our analysis of a similar Hamiltonian for our three-dimensional 
subsystem code more transparent. The Hamiltonian is given by nearest neighbor interactions constructed entirely 
from operators in the set T , 

n n— 1 

H = i Z i,j Z i,j+l + X 0,i X d+hi) ■ ( 27 ) 

z=l j=l 

Since H is constructed entirely from elements of T , this Hamiltonian can be decomposed as 

H = H s x tS z®I 2 . (28) 

s x ,s z 

To understand the exact nature of this Hamiltonian, we would need to diagonalize each of the H s x s z. What 
can be said, however, is that the ground state of the system will arise as the ground state of one or more of the 
H s x s z (numerical diagonalization of systems with a few qubits show that the ground state comes from only the 
sf = sf = +1 subsystem and we conjecture that this subspace always contains the ground state.) If k of the H s x s z 
contribute to the ground state, the degeneracy of the ground state will be 2k due to the subsystem corresponding 
to the C operators. Thus we see that we can encode quantum information into the subsystem degeneracy of this 
Hamiltonian, in the similar manner that information is encoded into the ground state of a Hamiltonian related to 
the toric code0, |3(| |33|. However, we do not know whether this system exhibits a gap in its excitation spectrum 
similar to that which exists in the toric codes. In Section IVII we return will introduce a similar Hamiltonian for our 
three-dimensional operator quantum error correcting subsystem. 



V. THREE-DIMENSIONAL OPERATOR QUANTUM ERROR CORRECTING SUBSYSTEM 

We now turn to a three-dimensional operator quantum error correcting subsystem which is a generalization of the 
two-dimensional subsystem code we presented above. In particular, whereas the construction for the two-dimensional 
model relied on the structure of T containing Pauli operators with even number of Pauli Z's in a row and even number 
of Pauli X's in a column, in the three-dimensional case we rely on a new set of operators with an even number of 
Pauli Z's in the yz plane and an even number of Pauli A's in the xy plane. 

Consider a cubic lattice of size n x n x n with qubits located at the vertices of this lattice and let n be odd. Let 
Oi,j,k denote the operator O acting on the qubit located at the (i,j, fc)th lattice site tensored with identity on all 
other qubits. We again use a compact notation to denote Pauli operator on our n 3 qubits by using two n 3 bit strings, 

n n ( Xij,k if Oi,j,k = 1 an d = 

P{a,b) = J] X%fZ°;f-*= [] \ Z ,< if Oij.k = and 6,^ = 1 , (29) 
i,j,k=i i,j,k=i { ~iYij : k if Oj^k = 1 and bij,k = 1 

where o,6 £ are n by n by n arrays of bits. 

As in the two-dimensional case, we will define three sets of operators, Tjj, S3, and £3 which are essential to 
understanding the subsystem structure of our qubits. 

The first set of Pauli operators which will concern us, T3, is made up of Pauli operators which have an even number 
of Xij_k operators in each xy-plane and an even number of Z^j^ operators in each yz-plane: 

{n n 1 

(-lfP(a, 6)|0 G Z 2j a i>jik = 0, and b iJtk = \ . (30) 
i,j=l 3,k=l ) 

These operators, like the analogous two-dimensional T form a group under multiplication. This group can be generated 
by nearest neighbor operators on our cubic lattice 

%, = \Xk,i,j Xk+isj, X^kjXi^+ij, Zij,kZij,k+i, Zi,k,jZi,k+i,j, Vi,j £ Z„,fc s Z n _i) (31) 
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The second set of Pauli operators we will be interested in is a subset of which we denote £3. 1S3 consists of 
Pauli operators which are made up of an even number of xy-planes made entirely of Pauli Z operators and an even 
number of yz-planes made entirely of Pauli X operators: 

{ n / n \ n I n \ j 

W)I0 A a ^,k =0,0 A =0 \. (32) 

S3 is an Abelian subgroup of T3 and all of the elements of £3 commute with all of the elements of 79. It can be 
generated by nearest xy-plane and yz-plane operators: 

/ n n \ 

= ( IT ^'U,fc^'M,fc+l, Y\. Zk,i,jZk+l,i,j>Vk S Z n _i \ . (33) 

We label these generators, as before: 

n n 

= Y\_ Xij^Xij^k, and Sf = Zij : kZi+ij.k- (34) 

i,j=l j,k=l 

S3 is again a stabilizer group. 

The final set of operators which we will consider, £3, is similar to £3 except that the evenness condition becomes 
and oddness condition 

!n j n \ n/n \ | 

(-l)*P(a,6)|0 e Z 2 ,0 /\ ai , jtk =1,0 A b ^ = 1 \ ■ ( 35 ) 
i=l \j.k=l J k=l \i,j=l / J 

£3 together with 1S3 forms a group and all of the elements of C3 commute with those of T3. 



A. Subsystem Structure 



All three of the sets, T3, S3, and £3 will play a directly analogous role to the sets T, S, and £ in our two- 
dimensional model. In particular if we let Ji denote the Hilbcrt space of our n 3 qubits, then we can partition this 
space into subspaces labelled by the 2(n — 1) different ±1 eigenvalues of the operators Sjf and Sf of Eq. |MJ Again 
we will label these eigenvalues by and sf , with s x and s z labelling these strings. The Hilbert space of the system 
then decomposes as 

H= H s x s z= H s x s z. (36) 

sf ,...,s*_ 1 ,sf ,...,s^_ 1 =±l s x ,s z 

Again, the 7i s x s z subspaces have a tensor product structure, 7i s x s z = T~tjx s z <8> T~tjx z, such that elements of £3 
act as 

L= I 2 „3_ 2n+1 ® L s x s z , \/L G £3, (37) 

s x ,s z 

and those of T3 act as 

T = T s x iSZ ®/ 2 ,vrer 3 . (38) 



B. Subsystem Quantum Error Correcting Procedure 



The subsystem error correcting procedure for the three-dimensional code nearly directly mimics that of the two- 
dimensional code. Here we discuss how the subsystem error correcting procedure works without going into the 
details as we did in the two-dimensional case. The three-dimensional procedure is nearly identical to that of the 
two-dimensional procedure with the sets T, 5, and £ interchanged with the sets 73,53, and £3 respectively. 
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We will again encode our quantum information into the Hf X — s z—i +1 \ n -i subsystem. Whereas for the two- 
dimensional code wc defined error strings for the rows and column conditions of the set T , now we define error 
strings for the xy and yz plane conditions of the set T%. If the Pauli error P(a, b) occurs on our system, then we can 
define the two error strings 

n n 

efc(a) = @ a<j,fe and /*(&) = £J) h,j,k- (39) 

i,j=l j-k=l 

Pauli errors with efc(a) = fi(b) = are errors from T3 and act trivially on the information encoded into the TL x _ z _ , . 

S k ~ S i — + 1 

subsystem. 

The quantum error correction procedure is then directly analogous to the one for the two-dimensional code. We 
measure the s£ and sf operators and treat these as nearest neighbor parity checks for a redundancy code on the 
fi(b) and efe(a) respectively. Then in direct analogy with the two-dimensional code, we can apply a subsystem error 
correcting procedure which restores the system modulo the subsystem structure. The three-dimensional code is a 
[n 3 , n, 1] code. We have thus gained nothing in terms of the distance of the code, but, as we will argue in the next 
Section, the three-dimensional code when converted to a Hamiltonian whose ground state is the subsystem code has 
intriguing features not found in the two-dimensional code. 

C. Logical Operators 

The logical operators for the three-dimensional code are directly analogous to those in the two-dimensional code. 
As in the two-dimensional code, operators from £3 can be used to enact Pauli operators on the information encoded 
into the T~ifx_ s z_ +1 subsystem. Similarly, we can enact a controlled-not between two encoded qubits by performing 

n 3 controlled-not gates between two identical copies of the code. This allows us to again perform any operation in 
the normalizer of the Pauli group on our encoded qubits. 

VI. SELF-CORRECTION IN THE THREE-DIMENSIONAL EXAMPLE 

In the 1930s, when Alan Turing wrote his now classic papers [38Ll39| laying out the foundations of computer science, 
there was absolutely no reason to believe that any computing device such as the one described by Turing could 
actually be built. One of the foremost problems, immediately apparent to the engineers of the day, was the lack 
of reliable components out of which a computer could be built. Von Neumann solved this problem, in theory, by 
showing that robust encoding of the classical information could be used to overcome errors and faulty components in 
a computer^?]]- Despite Von Neumann's theoretical ideas, it took the invention of the transistor and the integrated 
circuit, to mention only the broadest innovations, in order to bring forth the technological movement now known 
as the computer revolution. The overarching result of the technological innovations responsible for the computer 
revolution was the development of techniques which exhibited Von Neumann's theoretical ideas in a natural setting. 
Modern computers naturally correct errors in both the storage and manipulation of classical information. The task 
of robust storage and manipulation of the data is essentially guaranteed by the physics of these devices. There are 
distinct physical reasons why robust storage and manipulation of classical information is possible. 

If there are distinct physical reasons why robust storage and manipulation of classical information is possible, an 
obvious question to ask in the quantum information sciences is whether we can mimic these effects in the quantum 
domain. Do there exist, or can we engineer, physical systems whose physics ensures that the robust storage and 
manipulation of quantum information is possible? In this section, we will present evidence, in the form of a mean field 
argument, that a Hamiltonian related to the three-dimensional subsystem code might be exactly this type of system. 

A. Self-Correcting Quantum Memories 

The traditional approach to building a robust fault-tolerant quantum computer imagines building the computer 
using a complex microarchitecture of quantum error correcting fault-tolerant procedures. This poses a severe tech- 
nological overhead of controlling thousands of qubits in a complex manner, simply to get a single robust qubit. 
KitaevOli^ was the first to suggest that an alternative, less complex, method to constructing a fault-tolerant quan- 
tum computer might be possible. Kitaev showed that there exists a quantum error correcting code, the toric code, 
which is the degenerate ground state of a certain four body spatially local Hamiltonian on a two-dimensional lattice of 
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qubits. Kitaev imagined encoding quantum information into the ground state of this system and then, because there 
is an energy gap between the ground state of this system and the first excited state and because the errors which will 
destroy quantum information consist of error which scale like the size of the lattice, this quantum information would 
be protected from decoherence due to the environment as long as the temperature of the environment was sufficiently 
low. 

It is important to note that the Hamiltonian implementation of Kitaev's toric code (by which we mean encoding 
information into a physical system governed by the four-body Hamiltonian associated with the toric code), while 
providing a mechanism for the robust storage of quantum information, does not provide a full fault-tolerant method 
for quantum computation. The reason for this is that during the implementation of the physical processes which 
manipulate the information encoded into the ground state of the system, real excitations will be created which will 
disorder the system. This distinction has been confused in a manner because the toric codes can be used to construct 
a fault-tolerant quantum computer, but only with the aid of external quantum control which serves to identify and 
correct errors which occur during the manipulation of the information encoded into the system (see for example |36jp. 

The idea of a self-correcting quantum memory is to overcome the limitations of Kitaev's original model by con- 
structing a physical system whose energy levels not only correspond to a quantum error correcting code (in our case 
a subsystem code), but which also uses the energetics of this system to actively correct real errors created when the 
quantum information is being manipulated. Thus, while in the toric codes, a single real error on the system can 
create excitations which can disorder the system, in a self-correcting system, a single real error on the system cannot 
disorder the system. 

In order to explain the distinction of a self-correcting memory from the original toric code we will compare the 
situation to that of the one dimensional and two-dimensional classical ferromagnetic Ising model. These models will 
be analogous to the toric code Hamiltonian model and a self-correcting Hamiltonian model, respectively. 

Recall that in a ferromagnetic Ising model one takes a lattice of classical spins and these are coupled by Ising 
interactions between the neighboring spins via a Hamiltonian 



where Sj <E {±1} are the spin variables, the sum is over neighbors in the lattice, and J > 0. Notice that the 
ground state of this Hamiltonian corresponds to the uniform states Sj = +l,Vi or s, = — l,Vi. These, of course, are 
also the two codewords for a classical redundancy code. Thus we can imagine that we encode classical information 
into the ground state of this Hamiltonian in direct analogy to the way in which information (but quantum this time) 
is encoded into ground state of the toric code. Errors on the Ising codes are just bit flips. From hereon out when we 
refer to the Ising model we will implicitly be discussing the ferromagnetic Ising model. 

Recall some basic properties of the one and two-dimensional Ising models (see, for example flj ). We begin by 
discussing the thermal equilibrium values of the total magnetization, 



of these models. In one dimension, for any T > 0, the total magnetization of the Ising model vanishes in thermal 
equilibrium, whereas for the two-dimensional Ising model, the magnetization is zero above some critical temperature 
T c and below this temperature, two magnetizations of equal magnitude and opposite sign are maintained. Since the 
magnetization is a measure of the information recorded in the redundancy code, we see that if we encode information 
into the ground state of the one dimensional Ising model and this system is allowed to reach thermal equilibrium, then 
this information will be destroyed. On the other hand, for the two-dimensional Ising model, if we encode information 
into the ground state and the system is below the critical temperature T c , then this information will be maintained. 
Above T c , like the one-dimensional Ising model, the information will be destroyed. From the point of view of storing 
the information in the thermal states of these models, the two-dimensional Ising model is a robust medium, but the 
one-dimensional Ising model is not. 

But what about the properties of the Ising models on the way to reaching equilibrium (i.e. during the time evolution 
with the environment)? In the one dimensional case we find that the system will generically (depending on the exact 
method of relaxation) take a time which is suppressed like a Boltzman factor e~ J ^ T . Thus at low enough temperature, 
we can encode information into the ground state of the one dimensional Ising model and it will be protected for a long 
time. While the scaling of this decay rate is favorable in the temperature T, this type of approach is different from 
what is done in standard error correction where larger redundancy can be used to overcome errors without changing 
the error rate (as long as that error rate is below the threshold for the error correcting code.) 

What happens for the time evolution of a two-dimensional Ising model? If we start the system in one of the 
redundancy code states, then far below T c the system will relax quickly to the closest of the two equilibrium states 




(40) 




(41) 
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with a large total magnetization. As we raise the temperature closer to T c , this relaxation will slow down. Above 
T c the picture is similar to that of the one dimensional Ising model that if we are close to T c , then the relaxation to 
vanishing magnetization is suppressed like e~ J ^ T ~ Tc \ 

What are the main reasons for the differences in the ability to store information in the one and two-dimensional Ising 
models? A rough heuristic of what is happening is that in the two-dimensional Ising model, the errors self-correct 43J . 
Consider starting the one dimensional Ising model in the all Sj = +1 state. Now flip one of the spins at the end of the 
chain. This will cost an energy J. Flipping the neighbor of this spin will then cost no energy. Proceeding along the 
chain in this manner one sees that one can expend energy J to turn the system from the codeword all = +1 state 
to the all Si = — 1 state. Thus the environment need only supply this energy to disorder the system. However, in the 
two-dimensional Ising model, something different happens. Suppose again that we start in the all Sj = +1 state. Here 
if we flip a single spin (say on the boundary of the lattice) then the energy required to flip this spin is J times the 
number of bonds this spin has with its neighbors. Now flipping a neighbor will cost energy: in the two-dimensional 
model the energy cost of flipping a connected domain of spins is proportional to the perimeter of this domain. Since 
to get from the all Si — +1 codeword to the all Si = —1 codeword we need to build a domain of size the entire lattice, 
we see that we will require at least an energy times the size of the lattice to disorder the system. Suppose, now that 
errors are happening at some rate to all of the spins in the Ising models. In the one dimensional model, once one 
creates a single error, there is no energy barrier to disordering the system. In the two-dimensional model, however, 
there is now an energy barrier. In particular, the system coupled to its environment will not only cause errors, but 
will also cause errors to be corrected by shrinking the domains of flipped errors. As long as the error rate is not too 
strong (which corresponds loosely to being below the critical temperature T c ) the pathways that fix the error will 
dominate the actual creation of errors. 

Thus we see that a two-dimensional Ising model operating below the critical temperature is performing classical 
error correction on information stored in a redundancy code. In the one-dimensional Ising model and in the two- 
dimensional Ising model above the critical temperature, there is suppression due to a Boltzman factor, but there is 
no self-correction (or the self-correction is not fast enough) and the information stored in the redundancy code is 
destroyed. 

The two-dimensional anyon models of topological quantum computing and variations flol l3ol l37l El EL liH EiiL ItH 
S EES EHH including Kitaev's toric model, all share the property with the one dimensional Ising model 
that the system can disordered using only an energy proportional to the gap in the Hamiltonian. (The models in 
|50| contain errors similar to those in the two dimensional Ising model for certain types of quantum errors, but not 
for both phase and bit flip errors.) This can provide protection via an exponential suppression due to a Boltzman 
factor, but this does not provide indefinite correction. The idea of a self-correcting quantum memory, however is 
to mimic the two-dimensional Ising model. In particular below a critical temperature, quantum information stored 
in the system should persist even when the system is in thermal equilibrium and further the system should have a 
mechanism whereby real errors are corrected by the energetics of the system faster than the real errors occur when 
operating below the critical temperature. 

Finally, we note that there is one system which is widely suspected to be a self-correcting quantum memory. This 
is a version of the toric code on a four-dimensional lattice|36j. The problem with this model is that it exists in 
four dimensions and that it requires greater than two-qubit interactions in order to implement and is therefore not 
realistic for practical implementations. Our original motivation for considering the subsystem codes presented in this 
paper was to obtain a self-correcting system in the realistic setting of three or lower dimensions and using two-qubit 
interactions. 



B. The Three-dimensional Hamiltonian 



Next we turn our attention to a system which may be a self-correcting quantum memory. We provide evidence for 
this by showing that in a mean field approximation this Hamiltonian has properties for the expectation value of its 
energy which is similar to the energetics of the two-dimensional Ising model. 

Consider the following Hamiltonian on a cubic lattice of qubits constructed exclusively from elements of 7~s, 

n n—1 

H = —A [Xk,i,jXk+i,i,j + Xi^,jXi,k+ij + Zi y kjZi,k+ij + Zij^Zij^k+i) , (42) 

i,j=i fc=i 

with A > 0. This Hamiltonian consists of Ising couplings along the direction x in the xy-plane and along the direction 
z in the yz-plane. As in the two-dimensional Hamiltonian, Eq. I|27|l . we can use the fact that H is a sum of operators 
from 7j| to block diagonalize H with respect to the subsystem structure of our three-dimensional operator quantum 
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error correcting subsystem: 

H= ff,* a ,®Z a . (43) 

Information can then be encoded into the 1 2 subsystem. We would like to show that if we perform such an encoding, 
then the information stored in this subsystem will be protected from the effect of general quantum errors in a manner 
similar to that of the two-dimensional Ising model. 



C. The Mean Field Argument 

Ignore, for the moment the boundary conditions for the Hamiltonian and suppose that the ground state of Eq. I|42|l 
has the following properties for the expectation values of the Ising bonds in the system, 

(Xk,i,jXk+l,i : j)G = C X x 
{Xi y k,jXi^k+lj)G — c xy 
{Zi y k,jZ.^k+l,j)G = Czy 
{Zij,kZij,k+l)G = Czz, 

(44) 

where c a /3 > 0. In such a phase, the ground state energy is given by 

E G = (H) G = -Xn 2 (n - l)(c xx + c xy + c zy + c zz ). (45) 

Now consider the effect of a single Pauli error on a single qubit in the lattice (assume this is away from the boundary) . 
For example, consider a Pauli X error. Using the fact that X commutes with Ising bonds oriented along the x 
direction but anticommutes with Ising bonds oriented along the z direction, it then follows that expectation value of 
the energy of the system changes to 

Ei = (H) 1 =E G + 2\{c zy + c zy + c zz + c zz ), (46) 

where we have separated out each term arising from each of the four z direction Ising bonds which connect to the 
lattice site where the X error has occurred. Thus we see that a single bit flip error on the ground state will cause 
the expectation value of the energy to increase. Now consider the effect of applying a second Pauli X error which is 
a nearest neighbor to the original spin in the same yz plane. Now the Ising bond between these two errors does not 
contribute to the change in the expectation value of the ground state, but all of the z direction Ising bonds around the 
perimeter of the two flipped spins do contribute. Thus, for example, if the spins are neighbors along the y direction, 
the expectation value of the energy of this new state is 

E 2 = {H) 2 =E G + 2X(2c zy + Ac zz ). (47) 

Generalizing the above argument we see that a connected domain of Pauli X errors in an yz plane will result in an 
energy increase proportional to the perimeter of the domain. A similar argument will hold for Pauli Z errors in an 
xy plane, but now the Ising bonds along the x direction will contribute to the change in expectation value, and those 
along the z will not contribute. 

Suppose that a general error P(a, b) occurs on the ground state of H. Then in each plane yz plane the part of the 
error coming from a will produce excitations whose expectation of the energy scales like the perimeter of domains of 
errors in a and similarly for each xy plane, but now for the b component of the errors. From this argument we see 
that, at least for the expectation value of the energy, the system looks very similar to the two-dimensional Ising model. 
But now instead of only bit flip errors, more general quantum errors produce changes in energy of the system which 
are proportional to the perimeter of the erred domain. Thus we argue that this provides evidence that the model we 
have presented will be self-correcting. For the same reasons that the two-dimensional Ising model will not disorder 
the classical information stored in a redundancy code (i.e. since our errors require (expected) energy proportional to 
the perimeter of the erred domain) we expect the quantum information stored in our system will not disorder up to 
some critical temperature. 

Of course, the evidence we have provided is based on numerous assumptions arising from our mean field model. 
First of all we have ignored boundary conditions. It is possible that certain edge states could disorder the system. 
Secondly we have only made arguments about the expectation value of the energy after error have occurred to the 
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ground state. This doesn't give us concrete information about the energy level structure of our Hamiltonian. It could 
be that while the expectation value scales like the perimeter, there are actually error pathways whose energetics are 
much less favorable. Third we have assumed the existence of a phase with the desire expectation values of the bond 
energies. This phase may not exist, i.e. it may be that in the thermodynamic limit, the expectations values all vanish. 
This would totally invalidate the mean-field argument we have given above. Given these caveats our mean field 
argument only suggests that the system will be self correcting. Clearly rigorously establishing whether our memory 
is self-correcting is a challenging open problem. 

D. The Quantum Error Correcting Order Parameter 

In the Ising model, an indication that the information stored in a redundancy code is still there after thermalization 
is the persistence of the total magnetization of the system. In particular, for the two-dimensional Ising model at 
a temperature between zero and the critical temperature, the magnetization in equilibrium is never exactly equal 
to its maximal value, ±n 2 . This is because there are always small domains of flipped spins due to the interaction 
of the system with its environment. However, a magnetization which is different from zero may be interpreted as 
a measurement of the majority vote for the redundancy code in this system. If we are going to establish that our 
quantum system is actually self-correcting, it is important to identify an order parameter for our system which can 
be used to reveal the presence of quantum information in our system. This can be done using the operator quantum 
error correcting subsystem properties of the three-dimensional code. 

Suppose we encode quantum information into a quantum error correcting code and apply a number of quantum 
errors less than the number which the code has been designed to correct. We know from the theory of quantum error 
correction that the encoded information in this system can be recovered by the measurement of an appropriate error 
syndrome and the application of the appropriate recovery procedure. We can use this to construct an order parameter 
for any quantum error correcting code. 

A note about the nature of order parameters for quantum information before we describe this parameter for our 
three-dimensional model. In the classical Ising model, we found that below the critical temperature there was a 
bifurcation of the system into two magnetizations of equal and opposite magnitude. Since classical information is 
based upon a bit, we are not surprised to find that such a bifurcation into two states happens. Actually, depending on 
the initial state of the system before it is thermalized, the thermal state of the system can be any value between these 
two extremes. But if we start our system by encoding into one of the two codestates (the all ±1 state) then only the 
two bifurcated values will result after thermalization. What is the analogous situation for quantum information? For 
quantum information, we must show not that the bifurcation happens for a bit of encoded information but instead for 
a qubit of encoded information. Since a qubit is parameterized by the Bloch-sphere, one might expect that one needs 
an order parameter with similar properties. Such an order parameter can be constructed, but we can get away with 
examining fewer parameters in order to show the robustness of the quantum information. In particular if we make 
measurements along the x, y, and z directions of the quantum information, then because of the linearity of the density 
operator we can use this to show that the information has been preserved. In particular, we can imagine encoding into 
one of the eigenstates of Pauli operators along these directions and looking at this system after it has thermalized. 
Notice now that instead of a single order parameter, we will have three order parameters. In order to demonstrate 
the self-correcting nature of our three-dimensional Hamiltonian, we will need to show that the expectation value of 
these three order parameters each bifurcate below some critical temperature. 

Consider, now, the active recovery procedure for our three-dimensional subsystem code. We begin by measuring the 
and Sf operators. Given these measurements, we can, as in the active recovery procedure, deduce an appropriate 
recovery operator to restore the information originally encoded into a subsystem, modulo the subsystem structure of 
the system. If we were to apply this syndrome and measure the encoded Pauli logical operators for the code (which 
were given in Sec. IV C() then this would serve as an order parameter for our system. We note, however that if we are 
simply interested in measuring the Pauli logical operators and not in fully restoring the information into the original 
subsystem encoding, we do not need to actually apply the syndrome. This is because the syndrome we diagnose will 
either commute or anti-commute with the encoded Pauli operator we wish to measure. Thus given the syndrome we 
can, instead of applying the appropriate recovery operation, simply flip or not flip the answer we get from measuring 
the encoded Pauli operator depending on the exact syndrome measured. 

Now we can explicitly describe our order parameters. Suppose we measure the and Sf operators and obtain the 
values and sf and we measure one of the encoded Pauli operators. Suppose, for example, that this encoded Pauli 
operator is an encoded X operator which we measure by measuring all of the Pauli X operators in a fixed xy plane. 
Given the we can deduce whether an even or an odd number of Pauli Z errors will need to apply to the fixed xy 
plane to restore the quantum information in the subsystem (if possible). If this number is odd, then we simply flip 
the value of we measure for the encoded Pauli X operator and if this number is even, we do not flip the value for the 
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encoded X operator. Thus we see that we measure the x directional order parameter from our system by measuring 
the and the encoded X operator and, as a function of these values, produce a single number representing the x 
directional order parameter. Similar comments hold for the order parameters along the other cardinal directions. 

VII. CONCLUSION 

In this paper we have constructed a new class of quantum error correcting procedures based upon the notion 
of encoding quantum information into a subsystem. By encoding into a subsystem, we were able to demonstrate a 
recovery routine which explicitly used the subsystem structure. The three-dimensional code we constructed was shown 
to be related to a three-dimensional spin lattice system which we gave evidence for being a self-correcting quantum 
memory. We close by remarking on some open problems for this three-dimensional system and some thoughts about 
future directions for constructing self-correcting quantum systems. 

The first open question concerns the implementation of our model in a physical system. A particularly promis- 
ing system for such an implementation is with ultracold atoms trapped in an optical lattice[55|. Duan, Demler, and 
Lukin 56] showed how to simulate a large class of spin-spin interactions for these systems. An open question is whether 
their techniques allow one to implement our three-dimensional anisotropic spin-spin Hamiltonian. Of particular con- 
cern is the magnitude of the spin-spin coupling which one can achieve in these models. This will directly effect the 
critical temperature for any self-correction that occurs in the system. A further concern for the physical implementa- 
tion in an optical lattice is the ability to measure the syndrome operators S£ and Sf along with appropriate logical 
Pauli operators. Finally one would like to understand how to produce an effective controlled- not coupling between 
two such encoded lattices. Solutions to all of these problems would allow one to propose an experiment in which a 
self-correcting quantum memory could be demonstrated. 

A second open question is, of course, whether our three-dimensional system is indeed self-correcting. Noting that the 
Hamiltonian for this system does not possess a sign-problem, one approach to verifying this question would proceed by 
performing quantum Monte Carlo simulations of this system. The order parameters we have described could then be 
simulated at finite temperature and evidence for self-correction could then be examined. Another promising approach 
is to use recent new ideas in the density matrix renormalization groupie EH to simulate the thermal properties of 
this system. 

Another important question is whether one can design a self-correcting quantum systems in two dimensions. This 
would be particularly desirable if one wishes to physically implement the self-correction in a solid-state system. 

Finally a large open question is what role operator quantum error correcting subsystem codes can play in quantum 
information science. What do other such subsystem codes look like? Since these codes have important properties due 
to the fact that they exploit degenerate quantum codes, can subsystem codes be used to beat the quantum Hamming 
bound 32] ? Further, results which relied on showing that there were no subspace codes with certain properties (for 
example, as in [59^ need to be reexamined in light of the existence of operator quantum error correcting subsystem 
codes. 
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